Observation of Stiickelberg oscillations in dipole-dipole interactions 
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We have observed Stiickelberg oscillations in the dipole-dipole interaction between Rydberg atoms 
with an externally applied radio-frequency field. The oscillating RF field brings the interaction 
between cold Rydberg atoms in two separated volumes into resonance. We observe multi-photon 
transitions when varying the amplitude of the RF-field and the static electric field offset. The 
angular momentum states we use show a quadratic Stark shift, which leads to a fundamentally 
different behavior than linearly shifting states. Both cases are studied theoretically using the Floquet 
approach and are compared. The amplitude of the sidebands, related to the interaction strength, is 
given by the Bessel function in the linearly shifting case and by the generalized Bessel function in the 
quadratically shifting case. The oscillatory behavior of both functions corresponds to Stiickelberg 
oscillations, an interference effect described by the semi-classical Landau-Zener-Stiickelberg model. 
The measurements prove coherent dipole-dipole interaction during at least 0.6 ^s. 

PACS numbers: 32.80.Ee, 42.50. Hz, 32.80.Wr, 34.20.Cf 



I. INTRODUCTION 

The subject of the sinusoidal perturbation of a quan- 
tum mechanical object is a very general one and has a 
long history. Here we will treat this subject for the par- 
ticular case of particles with states that have a quadratic 
shift under the influence of a perturbation. This is what 
is generally found in second-order perturbation theory. 
The typical example is the polarizability of a neutral par- 
ticle under the influence of an electric field. The realiza- 
tion chosen to study this phenomenon is the dipole-dipole 
interaction between Rydberg atoms, because this system 
provides the resolution and flexibility required for these 
experiments. Apart from representing a much broader 
class of phenomena, the study of dipole-dipole interaction 
in oscillating fields provides us a direct tool to measure 
the coherence time of the interaction and the possibility 
of switching the interaction. Both aspects are relevant 
for the feasibility of implementing quantum information 
processing with Rydberg atoms. 

First we will give a description of polar and polarizable 
systems in a RF field, i.e. we compare the situation for 
states with a linear versus a quadratic Stark shift. We use 
the Floquet approach, which is applicable for periodically 
varying fields. In this framework a two-level system plus 
oscillating field is replaced by an infinite number of side- 
bands on the two states. For a quadratic Stark shift the 
result is fundamentally different and much richer than 
for a linearly Stark shifted state. Next we compare the 
Floquet description with the classical limit. In addition 
we interpret the oscillations of the population of the side- 
bands as Stiickelberg oscillations. 

In the subsequent section we describe details of the ex- 
periment and present results of the measurements. Re- 
lated experiments have been performed with transitions 
within single Rydberg atoms, in most cases for states 
with a linear Stark shift [U El |3] , but also quadratic shifts 
have been studied [1]. Collisions of Rydberg atoms in 



an oscillating field have been studied in an atomic beam 
setup [SJ 12] , as well as in a cold cloud [3 E] ■ In our exper- 
iment we control the interaction strength by maintaining 
a fixed distance between the Rydberg-atom volumes. Be- 
cause of the small unperturbed energy mismatch between 
initial and final states, we use a rather low frequency for 
the oscillating field; it lies in the radio frequency range 
rather than in the microwave range, giving the opportu- 
nity to test the coherence properties of the system for a 
longer period of time. 



II. A COUPLED TWO-LEVEL SYSTEM IN AN 
OSCILLATING FIELD 

To introduce the system under investigation we recall 
the two- level system [5J, coupled by an interaction V in 
the presence of a field F{t) 



F{t) ^ Fs + Frf cos Lot, 



(1) 



with an oscillating and a static part. We write the full 
Hamiltonian as 



H — Ho + Hp 



V, 



(2) 



where Hq is the Hamiltonian of the two-level system 
without interaction, the interaction with the static 
field and Hpup the interaction with the oscillating field. 
V is the dipole interaction between the two involved 
states. We want to know the population in both states, 
so we need to know the full wavefunction ip{r,t), by solv- 
ing the Schrodinger equation. We solve the problem by 
starting from the known static solutions for Hq + Hp^ , as 
described below, and then we add the other two terms as 



perturbation. In the next sub-section II A 'Numerical so- 
lutions in the time domain' we work out the problem by 
first adding the interaction term V to the Hamiltonian, 
and then the oscillating field Hp^^^. This provides us 
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with two coupled equations which can be solved numeri- 
cally in the time domain. In the subsequent sub-sections 
we work in the reverse order. In section [II B| 'Floquet ap- 
proach' we add the oscillating field Hp^^^ first and work 
out the problem in the frequency domain, which basi- 
cally provides us with an infinite number of steady states 
in the form of sidebands. The classical limit of the side- 
band population, section [TTC] gives some physical insight 
in its behavior. Then, in section II D 'Coupling between 
dressed and undressed state' we add the interaction term 
V to the Hamiltonian. The oscillations in the coupling as 
a function of RF amplitude can be understood in terms 



of Stiickelberg oscillations, explained in section HE 

We will start by treating the perturbation of the 
Hamiltonian in a static field. With just the first two 
terms of the Hamiltonian H = Hq + Hp^ we obtain the 
solution 



■0(r,i) = aijji{Y,t) + bip2{r,t), 



(3) 



where ipi and ^2 £^re the eigenstates of the system and a 
and b are their normalized time-independent amplitudes. 



V'i(2)(r,i) = ^i(2)(r)e-^'^i(^)*, 



(4) 



where Wi and W2 are the eigenenergies of the two sta- 
tionary states. 

We will discuss two cases, one for states with a lin- 
ear Stark shift and one for states with a quadratic Stark 
shift. The first case is depicted in Fig. [l^. State |1) al- 
ways has the energy and state |2) has the energy 
in the absence of any fields and a linear energy shift as 
a function of electric field F (all equations are in atomic 
units): 



W2^Wq- kF. 



(5) 



Here k is the permanent dipole moment of state |2). This 
system is in fact equivalent to two states with difference 
energy Wq and a difference in dipole moment of k, which 
makes this two-level system applicable to many realistic 
systems. The two levels cross at a field 



Wo 
k ■ 



(6) 



The second case is sketched in Fig.[T|3. Again, state |1) 
has the energy for all fields and state |2) has the energy 
Wo at zero field but the energy decreases quadratically 
as a function of electric field F: 



W2 = Wo 



-aF^ 



(7) 



Here a is the polarizability of state |2). Similar to the 
linear case this system is equivalent to two states with 
difference energy Wq and a difference in polarizability of 
a. The two levels cross at ii^quad, where 



F 



quad 



2Wo 



(8) 
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FIG. 1: (Color online): Energy level diagrams of state |1) and 
2) as a function of electric field. State |1) has no Stark shift, 
state 1 2) has a linear Stark shift in Fig. a) and a quadratic 
Stark shift in Fig. b). Due to the coupling between the states, 
the crossing becomes an avoided crossing, depicted in the in- 
set. When an oscillating field is applied, depicted as a dashed 
purple line, this translates as an oscillating energy for state 
1 2), depicted on the right hand side as a dash-dotted red line. 
The orange dotted lines depict the average energy during the 
oscillation and the corresponding field; in a) this field is sim- 
ply Fs, the average field, but in b) it is different and we define 
it as the efi^ective field _Feff, given in Eg. |25| 



A. Numerical solutions in the time domain 

The most direct method to calculate the evolution 
of the wavefunction in an oscillating field is to numeri- 
cally solve the two coupled differential equations for both 
states. For this we first add the coupling V and then the 
time dependent field perturbation Hp^^^ . When we add 
the interaction V , the degeneracy at the crossing is lifted 
and an avoided crossing arises, which has a width 
This r^o can also be interpreted as a quantum beat os- 
cillation frequency: the population of state V'l and V'2 
oscillate against each other with this frequency after an 
appropriate initial condition. The avoided crossing is de- 
picted in the inset of Fig. [T] where the width is given by 



r!o = 2(7/.2(r)|l/|V^i(r)). 



(9) 



When we add a time-dependent field it is convenient to 
separate the time and spatial dependence of the wave- 
function (Eq. [3]) 

V;(r,t) =Ti(t)^i(r)+r2(i)V'2(r). (10) 
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where '0i(r) and V'2(r) are the two original states with- 
out interaction. This separation of time- and spatial de- 
pendence is only allowed if both spatial states remain 
approximately the same despite the applied fields. We 



will see in section III A that this assumption is valid. Us- 
ing this wavefunction in the time-dependent Schrodinger 
equation we obtain 

*ri(2)(t) = M^l(2)(<)+ (11) 

Now the energies Wi and W2 depend on time. For the 
linear case we have 

W2 (t) = Wo-k{Fs + Frf cos ut) , (12) 

and for the quadratic case 

W2{t) = Wo - (F| + 2FsFrf cos ujt + F^p cos^ tut) 

(13) 

and in both cases Wi{t) = 0. 

We have numerically solved the two coupled equa- 
tions in Eq. |11[ However, the results are less informa- 
tive than the approach described in the next section, the 
Floquet approach. Although equivalent, an answer in 
terms of a coherent superposition of stationary states is 
more informative than the time-dependence of a wave- 
function. However, this numerical approach can be used 
for all time-dependent fields, whereas the Floquet ap- 
proach only works for periodic functions. 



B. Floquet approach 

In the Floquet approach we add the perturbations in 
the opposite order as in the numerical calculation; so we 
first perturb the states ipi and V'2 with the oscillating field 
and subsequently add the interaction V between the two 
states 111 [9]. As in the previous section we separate the 
time and spatial dependence, but this time we keep two 
separate wavefunctions 



V'i(2)(r,i) = Ti(2)(i) Vi(2)(r). 



(14) 



Using these in the timc-depcndcnt Schrodinger equation 
we obtain 

^f^2)it)^P^2)ir) = W^2){t)T^2){t)4>^2)ir). (15) 

We can factor out the spatial dependent part and inte- 
grate Eq. [15] 



(16) 



where Wi{t) = for both examples. 

For the linearly Stark shifted case, still following [9J, 
we use Eq. [12] to obtain 
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FIG. 2: (Color online)The sidebands of state |2) for the linear 
case in an oscillating field of 8 MHz for various RF amplitudes. 
The intensity of a line represents the population of the respec- 
tive sideband, where the original state, or the n = sideband, 
is colored blue and indicated with an " o" . The red dashed 
lines depict the borders of the classically allowed region. We 
used fc = 27r X 60 MHz/(V/cm) and Wo = 27r x 25 MHz. 



The last exponent can be written as an expansion of 
Bessel functions [10], such that 

^2(r,t)=V^2(r)e-^(^°-'=^«)* jJ'^U^^^' 



n— — 00 



This expression shows that the wavefunction consists of 
the original spatial wavefunction, but with a modified 
time dependence in the form of an infinite number of 
sidebands. This is the steady-state solution of the prob- 
lem, consisting of a number of states with different am- 
plitudes, given by the Bessel function J„ (^^^Jf^) , all with 
a different energy 



1^2 



Wo - kFs 



(19) 



g-i(VKo-feFs)tgi^i^ sinojt 



(17) 



In Fig. [2] the sidebands of ip2 are plotted for various RF 
amplitudes. It is clearly visible that more sidebands are 
populated when a larger RF amplitude is applied. The 
red dashed lines indicate the classical limits, the maxi- 
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FIG. 3: (Color online) The sidebands of state |2) for the 
quadratic case case in an oscillating field of 8 MHz for dif- 
ferent RF amplitudes. The intensity of a line symbolizes the 
population of that sideband. The n — sideband is depicted 
in blue and indicated with an " o" . The red dashed lines depict 
the position of the asymptotes of the distribution of occurring 
energies. We used the values corresponding to the experimen- 
tal values of the left resonance, as mentioned at Eq. |52[ 



mum and minimum field values that occur in the oscil- 
lation. This is equivalent with the absolute order and 
absolute argument of the Bessel function being equal. 

The result for the quadratically Stark shifted case is 
quite different. Combining Eqs. [16] and [13] we obtain 



T2{t) = e-^^o-s^C-FJ + s^^lF))* 



X e 



(20) 



Here we have two exponents to expand in terms of Bessel 
functions 



where n = m' + 2m is used on the third line. The am- 
plitude of the sideband, i.e. the inner summation in the 
above expression, is known as the generalized Bessel func- 
tion Jn{x,y). It was first used by Reiss [11] and investi- 
gated by e.g. Dattoli et al. p^. It is defined here as 



E 



Jn-2m{x)Jm{y)- 



(22) 



In practice the summation can be constrained to |m| < 
y + Sy^/"^ -I- 3 or |m| < \{n + x + 3x^/^ + 3), because be- 
yond these points the Bessel function is negligibly small 
[I3I . Instead of evaluating the sum, it is also possible to 
calculate the generalized Bessel function on basis of the 
recurrence relation [14 in analogy with the well-known 
approach to calculate the regular Bessel function on basis 
of its recurrence relation. 

Using these results the wavefunction becomes 



( aFRpFs aFl^p 



00 
T, — — 00 



* 7 



8uj 



(23) 



A similar description of sidebands on quadratically 
shifted states is given in |3]. Note that the leading order 
of the series expansion is proportional to F^p] a result 
that is also found for the linear case and that is expected 
from a description of multi-photon transitions. 

If we compare this wavefunction with the linear case it 
is the generalized Bessel function instead of the regular 
Bessel function that determines the sideband amplitude. 
The result depends on both the static field Fs as well 
as the RF amplitude Fj^p, whereas for the linear case it 
depends on the RF amplitude only. The energy of the 
sidebands is also fundamentally different and one obtains 



W2^n = Wo- -a 



RF 



(24) 



Whereas in the linear case the energy of a particular side- 
band shifts with the static field Fs only, in the quadratic 
case there is an additional term that depends on the RF 
amplitude \ctF^p. This term corresponds to the AC- 
Stark shift of the state. It is convenient to define an ef- 
fective field Foff , such that the energy shift of such a static 
field corresponds to the energy shift of a given static plus 
RF field: 



Fi 
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:F] 



RF- 



(25) 
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Y Jn-2mix)Jm{y) , (21) 



\m— — oo 



The different behavior is visible in Fig.[3| where the side- 
bands are plotted versus static field for several RF am- 
plitudes. Again, the intensity depicts the population of 
the sideband: Jn {aF]ipFs/uj,aF^p/{8uj)) . As in the 
linear case more sidebands are populated with a larger 
RF field, but a clear difference is that the population of a 
sideband is now not only a function of the RF amplitude. 
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but also of the DC field, visible as the changing intensity 
within one sideband. Secondly, the difl^erent positions of 
the sidebands in each subfigure illustrates that the side- 
band energy also depends on Fj^p (note that the original 
state, or the n = sideband, in blue and indicated with 
" o" , is shifted downwards for larger Fnp ) . 

The next step would be to add the coupling between 
state |1) and state |2) (section II D I, but first we will look 



closer into the population of the sidebands for different 
frequencies. 



Classical limit of the sideband population 
distribution 



so there are not many populated sidebands. For smaller 
frequencies (or large x,y), more and more sidebands arise. 
Around n — the population oscillates and near \n\ = \x\ 
the Bessel function reaches a maximum, after which it 
drops exponentially to zero. Interestingly, the sideband 
populations extend further in energy beyond \n\ > for 
the higher frequencies as for the lower frequencies. For 
the number of populated sidebands beyond \n\ = \x\ the 
opposite is true. These extensions can be understood as a 
quantum mechanical effect, whereas the energy levels are 
in a classically forbidden region. The generalized Bessel 
function also has these features. In addition, it has a 
more asymmetric structure, which we will discuss further 
below. 
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FIG. 4; (Color online) Depicted in black is the sideband pop- 
ulation of state j2) as a function of the sideband energy (Eqs. 
19 and 241 for various frequencies, u is 2n times the denoted 
frequency. We used again fc = 27r x 60 MHz/(V/cm) and 
Wo ^ 2tv X 25 MHz and field values of Fs = 0.2 V/cm and 
Frp — 0.3 V/cm. The red dashed lines depict the distribution 
of clasically occurring energies, calculated with the approach 
of Eq. |26] 

To illustrate the sideband population and to obtain 
some physical insight into its behavior, we depicted the 
sideband population for a fixed field amplitude and offset 



as a function of the sideband energy (Eqs. 19 and 24) 
and various frequencies in Figs. |4] and |5] So we plot the 
squared (generalized) Bessel function J„(x)^ respectively 
Jn{x, yY for fixed x and y (per subfigure) as a function of 
n. The pictures can also be seen as a vertical cut through 
a subfigure of Fig.|2]or[3]with the appropriate parameters 
Fg, Ffip and uj. For large frequencies, i.e. small x and y, 
the function peaks at n = and drops off on both sides. 
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FIG. 5: (Color online) Similar as Fig. |4] We used the experi- 
mental values for the left resonance (Eq. |52[ ) and field values 
of Fs =0.2 V/cm and Frf =0.45 V/cm. The distribution of 
occurring energies (red dashed lines) is calculated with Eq. 26 
Note the occurence of three turning points. 

For low frequencies the (generalized) Bessel function 
approaches the classical limit, meaning it approaches the 
density distribution of the occurring energies. To calcu- 
late such a density distribution it is useful to know how 
a distribution, /(x), changes under a change of variable, 
from x to y. li y — g{x) the resulting density function 
becomes 



p{v) - E 



1 



9'{9l\v)) 



(26) 



with g' the derivative and g^^ the inverse function and 
n the number of solutions in x for g{x) — y. g^^{y) are 
these solutions. 

We could choose g to be the energy as a function of 



time, given by Eq. 12 and 13 It is, however, instructive 
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FIG. 6: (Color online) As Fig.jsj but at a lower frequency and 
zoomed in around the middle asymptote. The blue dotted 
curve depicts the moving average of the sideband population 
depicted in black. 



and easier to first calculate the distribution of the occur- 
ring fields. So we use g{t) — F{t) — Fs + Fupcosajt. 
We start with / as a uniform normalized distribution. It 
is sufficient to have instances during half a cycle of the 
cosine, 



f{t) = 



^ if < Cji < TT 

7T — — 

elsewhere 



(27) 



Apphcation of the functions / and g to Eq. [26] provides 
the probability density distribution of the field 



1 







if Fs-Frf<F<Fs+Frf 
elsewhere 



(28) 



This function diverges at F = Fs±Frp. The asymptotes 
correspond to the turning points of the field. 

The function F{t) is sketched in the energy level di- 
agram Fig. [T] as a dashed purple line. How the field 
oscillation translates to an oscillation in energy (Eq. [T2| 
and 13) is depicted on the right hand side as the red 
dash-dotted line. In the linear case, the function sim- 
ply remains a cosine, but in the quadratic case the func- 
tion is distorted. For Frf > Fs the density distribution 
consists of two parts that add up. In this case there is 
a third turning point when the field goes through zero. 
Both functions are plotted together with the strength of 
the sidebands in Figs. |4] and [S] Here they are multiplied 
by the frequency to match the average height of the side- 
bands. The asymptotes of the distribution functions are 
also depicted in Figs. |2] and [3] as red dashed lines. 

We observe in Figs. |4] and 5 that for the lower frequen- 
cies the sideband population indeed converges nicely to 
the distribution functions. The height of the peaks os- 
cillates around the red dashed line. In most cases values 
occur between and twice the red dashed line. In the 
case of the generalized Bessel function, on the right hand 
part of the plot some much higher peaks occur. However, 
the moving average still approaches the red dashed line. 
This is depicted in Fig. [6] where we zoomed in around 
the asymptote Wq — ^a{Fs — FrfY- 
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FIG. 7: (Color online) Calculated interaction strength as a 
function of static and oscillating field for the linear case. The 
vertical lines depict the nuj resonances, where n varies from 
—6 (the one most right) to 6 (the one most left) and u is 
fixed to 27rx4 MHz. The darker the color the stronger the 
interaction. Destructive interference between the two states 
is depicted by the curved lines in red and green, according to 
Stiickelberg theory described in subsection |II E| The dashed 
blue lines give the boundaries of the classically allowed region; 
in other words these lines depict field combinations where the 
crossing (Eq. |6| occurs at an extremum of the RF oscillation. 

Now that we know the wavefunctions in an oscillating 
field, we can add the coupling V between tpi and "02 ■ Ev- 
ery sideband of '02 crosses with ipi . After introducing the 
coupling V , every crossing becomes an avoided crossing. 
When we use the wavefunctions given by Eq. [14] we get 

(7^1 (r, t)\V |02(r, t)) = (7^1 (r)| V T^T^. (29) 



The n*"^ sideband of the linear case (see Eq. 18 1 becomes 
(Vi(r,i)|T/|V^2(r,t)), 



(0i(r)|y|02(r)) 
fkFRF_\ 



The resonance condition is seen in the exponent 
nw = Wo — kFs- 



(30) 



(31) 



The coupling strength at this resonance (the half width 
of the n**^ avoided crossing) is given by the time indepen- 
dent part of Eq. [30] 



nr. 



= (^i(r)|y|^2(r)) J„ 
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(32) 
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FIG. 8: (Color online) Calculated interaction strength as a 
function of static and oscillating field for the quadratic case. 
The quarter circles in varying gray depict the no; resonances. 
Here n runs from —25 (in the upper right corner) to 3 (in 
the lower left corner); the frequency of the RF field is here 
oj = 27rx8 MHz. The darker the color the stronger the inter- 
action. Destructive Stiickelberg interference between the two 
states is depicted by the curved lines in red, green and orange 
(subsection |II E[ ) . The dashed blue lines give all field combi- 
nations where the crossing is reached during an extremum of 
the RF oscillation. 



Here we used Eq. |9] which gives the coupling strength 
fio without oscillating field. So the total coupling ^fio 
is distributed over all avoided crossings, according to the 
Bessel function. In a similar manner, we get the reso- 
nance condition for the quadratic case 



nuj ^ Wo - -a \ Fg 



1 



RF 



and the coupling strength is 



~2' 



_ ( aFRpFs aFlp 



(33) 



(34) 



In Fig. [7] we plot the coupling resonances as a function 
of RF amplitude, Fjip, and static field, Fg, for a fixed fre- 
quency uj for the linear case. The grayscale depicts the 
squared coupling strength or Jn (kFjip/uj)'^ . The reso- 
nances are just vertical lines, since the resonance depends 
on the static field only. The blue dashed lines indicate 
the boundaries of the classically allowed region, which 
correspond to the regions where the field F{t) (Eq. [T]) 
reaches the static field resonance (Eq.[6]). These clas- 
sical boundaries correspond to the boundaries depicted 
as dashed red lines in Fig. [2] and to the asymptotes de- 
picted as red dashed lines in Fig. [4j Above these lines the 



coupling strength oscillates as a function of Fjip, which 
we can interpret as Stiickelberg oscillations, described in 
the next section. 

Fig. [8] shows the corresponding plot for the quadratic 
case. Here the resonance lines are curved, in fact, each 
line follows a line of constant effective field as defined 
in Eq. 25 The boundaries of the classically allowed re- 



gion are again indicated with the two lower blue dashed 
lines. They correspond to the lines Fs+Fjip = Fq^g^^ and 
Fs-FjiF = Fquad (with Fquad defined in Eq. [s]). There is 
a third blue dashed line at higher RF amplitudes, which 
indicates the boundary Fs — Frp = —Fquad; above this 
line both the crossing at —Fquad as well as the crossing 
at +Fquad are involved. Note that the third boundary 
or asymptote from Figs. |3] and [5] the one with value Wq, 
is not relevant, because it is never involved in the cross- 
ing with state |1). Within the classically allowed region, 
i.e. above the lower two blue dashed lines, we see the 
Stiickelberg oscillations. 

At this point there is enough information to start the 
measurements; it is known where to expect resonances 
and how strong they are. However, we will first explain 
why the oscillations in the interaction strength can be 
interpreted as Stiickelberg oscillations. 



E. Stiickelberg oscillations 

Stiickelberg oscillations [T5] are due to an interference 
effect between two quantum states, which occur when 
an avoided crossing is traversed at least twice. Imag- 
ine we start with all population in the state |1) and 
subsequently we sweep the field through the resonance 
and back (Fig. [T|). When sweeping through, part of the 
wavefunction will follow the state |1) diabatically and 
'ignores' the avoided crossing, while the other part will 
follow the eigen-energy adiabatically and end up in state 
1 2). When coming back again to the crossing, both parts 
of the wavefunction have built up a different phase, due 
to their energy difference and due to a difference in pass- 
ing the crossing. This phase difference leads to interfer- 
ence, which can be observed when varying for example 
the duration or the amplitude of the field sweep. In this 
particular case we have a sinusoidal oscillation of the field 
with a fixed frequency and we vary both Fg and F^p. 

To calculate the interference patterns we use the 
Landau-Zener-Stiickelberg (LZS) model [HI [TBI El) as 
described in e.g. [TB] and apply it to the situation 
sketched in Fig.[T| It is assumed that the avoided crossing 
(Eq. |9]) is small compared to the range of occurring ener- 



gies during an RF oscillation (Eq. 12 or 13). In that case 



we can separate the acquired relative phase into a phase 
evolution due to energy differences and phase jumps due 
to avoided crossings. The Hamiltonian of the two-level 
system in a diabatic description is 



H 



Wait) V 

V* Wb{t) 



(35) 
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where Wa and Wt are the unperturbed eigen-energies 
of respectively state \a) and \b), and V is the coupling 
between them. In order to better exploit the symmetry 
of the problem, we choose a definition of energies that is 
slightly different from before: Here Wb — —Wa = ^^^2, 
with W2 defined in Eq. [l2] or [13] The energies Wa and 
Wt during an RF oscillation are illustrated in Figs. [9] and 
[TT] for a specific choice of Fs and Fj^p . 

In the LZS model we describe the evolution of the 
wavefunction with transfer matrices. The transition at 
a level crossing is described by the transfer matrix (TH] 



M = 



i(f) 



(36) 



The matrix M is applicable for crossings where state \a) 
crosses \b) from the lower-energy side. In the opposite 
case the transpose of M, should be used. The phase 
angle (j) is the so-called Stokes phase [T^ 
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FIG. 9: (Color online) The energies of states |a) and \h) as a 
function of time during the RF-oscillation for the case of linear 
shifting Stark states. Indicated are the times where the two 
states cross, ti, and ta, and the matrices that describe the 
evolution of the wavefunction at the indicated time intervals: 
M^, G12, M and G23. The green and red colored areas 
indicate the phase evolution of the state which is 

contained in the Gi, matrix. 



with values between and 7r/4. And e is the Landau- 
Zener adiabatic transfer probability 



1 - e 



(38) 



Both quantities depend on the adiabaticity parameter b 
which in turn depends on the coupling strength V and the 
rate of change of the difference of the diabatic energies 



\d{Wb-Wa)\ldf 



(39) 



In the diabatic limit & is small and subsequently e is small 
as well, while (/) approaches the value 7r/4. In the adia- 
batic limit e approaches 1 and (/> approaches 0. 

After the crossing, the fraction in state |6) develops a 
different phase compared to the fraction in state |a), be- 
fore coming back to the crossing. Both phase evolutions 
are described by the matrix 







(40) 



with = ^Wb{tY ^V^dtK Wb{t)dt. The ap- 
proximation is valid for small V . Note that is a 
function of Fg and Frf- The matrices M and G which 
describe the phase development at and between the cross- 
ings are indicated in Figs. [9]and[TT| The shaded surfaces 
indicate the phase of state |5) as it develops between the 
crossings. 

For states with a linear Stark shift the situation is de- 
picted in Fig. [9] and, assuming that wc start with all 
probability in |a), the amplitude of both states after N 
oscillations is given by 



T\N 



(GzsMGiaM^) 



= S 



N 



(41) 



where S describes the time evolution of the wavefunc- 
tion over one RF oscillation. The matrix S is unimodu- 
lar (det S — 1). Therefore we can write S'^ in terms of 
Chebychev polynomials J7„(C) [20] 



S^^ =S[/Ar_i(e)-I(7jv-2(C), 



(42) 



where f = Tr(S). The solution of the full wavefunction 
in the linearly Stark-shifting case after N oscillations can 
be expressed in terms of the above defined quantities. 
In view of the size of the expression, we just give the 
population of state |6) 



Pfc = 4(l-e)esin2(ei2 + 0) C/^_i(0 



with the trace 



i = cos(ei2 + (/)) cos(e23-(?!') 
- (l-2e) sin(ei2 + 0) sin(e23^ 



(43) 



(44) 



In the near diabatic limit (e<C 1) the population Pf, re- 
duces to 



Pb = 4esin2(ei2 + ?!)) 



sin'(jV(ei2 + e23)) 

sin2(ei2+e23) 



(45) 



In Fig. 10 we depicted the population of state |&) as a 
function of the static field Fs and the RF amplitude Frj? 
for = 3. Note the similarity with Fig. [7] For larger A 
the vertically elongated areas of large population become 
narrower. 

The quadratic case allows for a similar approach, but 
there is a twist. In cases where the states cross each other 
at only one field, either -l-Fquad or — Fquad (Eq- [8|, the 
same equations are used as for the linear case. However, 
for large RF amplitudes it is possible that the two states 
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FIG. 10: (Color online) The population of state \b) as a func- 
tion of the static field Fs and the RF amplitude Frf (Eq. 
45 I for A'' = 3, where red is a large population and blue is a 
population of 0. 
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FIG. 11: (Color onhne) Similar to Fig. [9] now for the case of 
quadratically shifting Stark states and large RF amplitudes, 
such that both the -|--Fquad and the —Fquad crossing occur. At 
ti and t4 the field is at the -|-Fquad crossing and at t2 and ta 
the field is at the — -Fquad crossing. Furthermore, it is clearly 
visible that 0i2 and 634 are equal, they cover the section 
between +i^quad and — i^quad- 



cross at two diflFerent fields: F — ±_Fquad- The speed 
at which each crossing is crossed is usually different, so 
we have different values for the adiabaticity parameter 
6 which will be labelled + and — . The same holds for 
parameters that depend on <5: 0, e and M. The situ- 
After N oscillations the 



ation is illustrated in Fig. 11 

wavefunction has evolved to 



T\N 



(G45M^ G34M^ G23M_ Gi2M^ ) 



(46) 



The parameter of the Chebychev polynomial now be- 
comes 

e = (i-ejcos(2ei2+e23+e45) 

+ e_ cos(2(t)_ +2612-623+645) 

+ 2sin(645-0+) e+ sin(26i2+623 + <?!>+) 

- 2sin(623-0_)(e+e_ cos(26i2 + 0_ +0+) 

(47) 



+ J(l-e_)e_(l-eje^ 



where we have used 612 = 634. After N oscillations the 
1 6) population is 



Pb =4 (l-2ej./(l-eje_ sin(623-0J 



(l-e+)e+(sin(26i2+623 + 0+) 

-2e_cos(26i2 + 0++0jsin(623-0j)]V^_i(O 

(48) 

For near-diabatic transitions e_ , ^1 the transferred 
population reduces to 



Ph =4 



sin^(jV(26i2 + 623 + 645)) 

sin2(26i2 + 623 + 645) 



i7sin(26i2+623 + <^+)-y^sin(623-(/'J • (49) 



This \b) population is plotted as a function of the static 
field Fs and the RF amplitude Frf for TV = 3 in Fig. [12] 
This figure shows a similar structure as Fig. [8] 

In the expressions for the probability of state \b) (Eqs. 
[45] and 49 1 we recognize two factors. Firstly, the familiar 
shi^{NxJ/ sin^(a;) factor which gives rise to the resonance 
condition x = riTr; a very ubiquitous phenomenon. Work- 
ing out this resonance condition we get exactly the same 
result as before: Eqs. |3l] and |33) Secondly, we identify 
what we call the Stiickelberg interference factor: 



Ii = 4esin^(6i2 



(50) 



for the linear case or the quadratic case where only one 
crossing occurs and 



h = 4[^e_^ sin(26i2+623 + <^+) - y/FZsin{e 



'23- 



(51) 

for two-crossings case. Note that these interference fac- 
tors are independent of N and that they occur already 
for one single oscillation. 

In Fig. [7] we have depicted the regions of destructive 
interference. The green lines are the destructive inter- 
ference regions for the situation described above. The 
red lines give the destructive interference regions for the 
case of a different starting point of the RF oscillation: 
between ti and ^2 (Fig. |9]). Note that in this case the 
phase difference that results in the interference (Eq. 50 1 
involves 623 (the red area in Fig. |9| instead of 612 (^the 



10 




0.2 0.4 0.6 

Static field (V/cm) 



0.8 



FIG. 12: (Color online) The population of state \b} as a func- 
tion of the static field Fs and the RF amplitude Frf (Eq. 
49 1 for A'' = 3, where red is a large population and blue is a 
population of 0. 



green area). This is because only the phase difference 
that has developed between the first two crossings is rel- 
evant; the last part gives just an overall phase factor to 
the wave function. 

In Fig. |8] the double crossings case occurs in the upper 
left corner, above the upper blue dashed line. Below this 
blue dashed line (but above the lower two blue dashed 
lines) the one-crossing case is valid. In this lower part we 
have equivalently plotted the destructive interference in 
red and green. For the two-crossings case the interference 
factor I2 (Eq. 50) involves the first three colored area's 
O23 and 034(which is equal to O12). 



(Eq. 

from Fig. [Til 



612 



So here the phase difference that has developed between 
ti and ^4 is relevant. The orange regions in Fig. |8] corre- 
spond to the regions where this factor /2 = 0. Starting 
between t2 and ts, the interference factor I2 would con- 
tain instead of 023 and the relevant phase difference 
has developed between ^3 and tg (ig is the moment of 
the crossing after t^); this is indicated as the red regions 
in Fig. |8] Other starting points are possible, but not 
indicated in the figure since no new physics is involved. 

The semi-classical Landau-Zener-Stiickelberg (LZS) 
model we described is an alternative to the Floquet 
model, described in sections II B and II D[ Both mod- 
els give approximately the same results, as we compare 
Fig. [To] with Fig. [7] and Fig. [12] with Fig. [5] Secondly, 
the regions of destructive interference correspond nicely 
with the minima of the (generalized) Bessel function in 
Figs. [7] and ]8] This shows that both theories work well 
in the description of our measurements. The Floquet ap- 
proach is valid for an infinite number of RF periods. In 



the experiment we work with an 8 MHz oscillation dur- 
ing 20 /iS, i.e. 160 periods, so we believe that the Floquet 
approach is very well validated. The LZS description is 
not fully quantum mechanical - the exact time evolution 
near the crossings is not taken into account - and it does 
not give an answer outside the classically allowed regions, 
i.e. below the lower blue dashed lines in Figs. ]7] and ]8] 
However, the LZS description might give more physical 
insight into the problem, such as the interpretation of 
the occurrence of maxima and minima in the interaction 
strength. 



III. EXPERIMENT 

A. Spectroscopy 

In the following part we will describe the performed 
experiments and compare the results with the theory pre- 
sented above. For the two interacting states we have used 
two two-atom states, which are coupled through dipole- 
dipole interaction. We prepare atoms in the 41d state 
in a confined region of space, and some tens of microns 
away, we prepare atoms in the 49s state; in both cases 
by laser excitation. At a static field of approximately 
0.4 V/cm the binding energy of one 41d3/2 atom plus 
the binding energy of one 49si/2 atom equals the sum 
of binding energies of one 42pi/2 atom and one 49p3/2 
atom. The two-atom energy level diagram is depicted 
Both the 41d3/2 and the 49p3/2 atoms un- 
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in Fig. 

dergo Stark sphtting of the two \mj\ states, giving four 
two-atom states. The 42pi/2,i/2 + 49p3/2,i/2 state (no- 
tation: ntj\m^\) is plotted as the dotted red line and 



the 42p 



1/2,1/2 



49p, 



3/2,3/2 



state is plotted as the dash- 



dotted green line. The 41d 



3/2,1/2 



49si 



/2,l/2 



state is 



shown as the solid blue line. The 41d3/2,3/2 state is not 
excited by the laser, so the 41d3/2,3/2 + 49si/2.i/2 state, 
depicted with a blue dashed line, does not play a role 
in our experiment. The energies and polarizabilities are 
calculated with the Numerov method ^1] and we get 
Wo =25.15(13) MHz. The transition 

41d3/2,l/2 + 49Si/2,l/2 ^ 42pi/24/2 + 49p3/2,l/2 (52) 

has a difference polarizability of ai — 347.04(4) 
MHz/ (V/cm) ^ and is therefore resonant at a field of Fi 
= 0.3807(15) V/cm. The transition 

41d3/2,l/2 + 49Si/24/2 <-> 42pi/24/2 + 49p3/2,3/2 (53) 

has a difference polarizability of a2 — 297.40(4) 
MHz/ (V/cm) ^ and is resonant at a field of F2 = 
0.4113(16) V/cm. 

The states couple through resonant dipole-dipole in- 
teraction. The single-atom states have virtually no per- 
manent dipole moment, but one atom oscillates due to 
the coherent 49s and 49p superposition and the other 
atom due to the 41d and 42p superposition. These su- 
perpositions of two states do have a dipole moment, a so 
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FIG. 13: (Color online) The two-atom energy levels of the 
41d-f 49s and 42p-|-49p system. In the inset we have zoomed 
in around the two relevant crossings, which become avoided 
crossings due to the dipole-dipole interaction. 



called transition dipole moment /x, fi^^ — (49s|r|49p) and 
t^dp = (41d|r|42p). Transition dipole moments oscillate 
with a frequency that is given by the energy difference of 
the contributing states. Both dipole moments are around 
lOOOooe but can have various polarizations: fj, — for 
Amj = or = /Za; ± i^y for Arrij — ±1. Note that, 
in contrast to the situation in a constant field, here the 
sign of rrij is relevant. The strength of the dipole-dipole 
interaction is given by the well-known expression 

Vdd = (54) 

and is approximately equal to 2tt 100 kHz for a distance 
of 25 fim. This results in typical values of the adiabaticity 
parameter 6 between 10~^ and 10~^ (depending on the 
values of Fs and Fjip). This means that the crossings 
are traversed in a highly diabatic manner. The adiabatic 
transfer probability e is between 10^"* and 10^^, so the 
approximations eqs. |45] and |49] are well justified. The 
Stokes phase ip lies between 0.2487r and 0.257r. 

All 5 involved Rydberg states have a purely quadratic 
downwards Stark shift. This energy shift is due to the 
coupling with several £ states in a field. This means for 
the eigenstates in a field £ is not a good quantum number 
and the orbital angular momentum is not a conserved 
quantity. It is then not allowed to split the spatial and 
temporal dependence of the wavefunction (as in Eq. 14 1, 
because the field changes with time, and therefore 
changes with time. This is not a problem with linearly 
Stark shifted states, because here these states have fc as a 
good quantum number, the parabolic quantum number 
which determines the dipole moment of the state. So 
in the linear case iIj{v) does not change with the field 
(and therefore not with time). Fortunately, we need to 
use only very small fields. In the field of 0.4 V/cm we 
calculated the eigenstates of a whole range of Rydberg 
states and it appears that all involved states are more 
than 99.8% pure, e.g. the population of the 49p part of 
the wavefunction we call 49p is 99.8%. With this number 



we are confident that the separation of the space and time 
dependent part of the wavefunction is justified (see Eq. 

To|. 

In the experiment we explored the coupling strength 
as a function of the applied static field and the intensity 
of the radio frequency wave for this quadratically shifted 
system. In fact we made a few cuts though Fig. |8] a few 
horizontal scans, one scan along the vertical axis and fi- 
nally several scans, where we followed a curved resonance 
line and could clearly observe the oscillations. The results 
(section HID) are discussed after we have described the 



experimental setup in section |IIIB| and the simulations 



we have done on a system of multiple atom pairs in III C 



B. Experimental Setup 

To control the dipole-dipole interaction strength Vdd 
(Eq. 54) we confined the atoms in space by using focused 
laser beams for the creation of the Rydberg atoms, as 
described in Two sets of pulsed dye-lasers excite 

ground-state atoms in a magneto-optical trap (MOT) to 
the Rydberg states. By using cold atoms, we ensure that 
the atoms hardly move on the timescale of the experi- 
ment. The cylinder-like volumes that contain the Ryd- 
berg atoms are approximately 15 iim in diameter, deter- 
mined by the laser focus. The length is about 0.5 mm, 
determined by the diameter of the MOT cloud. Each vol- 
ume contains about 25 interacting atoms. The distance 
between the two laser beams is chosen to be 25 yum. Sev- 
eral ms before the laser excitation the magnetic field for 
the MOT is switched off, because this field broadens the 
resonances due to position-dependent and nij dependent 
variation of the Zeeman effect. 

The atoms are located between two field plates, on 
which the 8-MHz oscillating field is applied as well as an 
ionizing field ramp to detect the atoms and their states 
[9]. The oscillating field is applied by an Agilent 33250A 
Arbitrary Waveform Generator, controlled by an NI Lab- 
view program. The distance vector between the cylinders 
is in the same direction as the electric field. We use the 
fraction of 49p atoms {N4gp/ [A^49s -I- iV49p]) as a measure 
for the interaction strength. The exact field is calibrated 
with a static field scan over the known resonances Fi 
and i^2- The atoms are detected after 20 /iS of interac- 
tion time. All data are averaged over 200 realizations. 
More details on the experimental setup can be found in 
[231, m and I221- 



C. Simulations for multiple atom pairs 



In section [IIP I we derived expressions for the coupling 
strength as a function of the oscillating field. However, 
what we measure in the experiment is the 49p fraction, 
which is a measure for the 49p-|-42p state (pp-state) pop- 
ulation. The population of this pp-state on resonance 
is given by Ppp = sin^rit/2, since we start with 100% 
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in the sd-state. The interaction strength fl is given by 
Eq. 34 which contains the interaction V (Eq. |9|. This 



interaction is in this case the dipole-dipole interaction 
strength Vdd (Eq. 54 1, which depends strongly on the 
atom- atom distance, which is restricted but not fixed. 
To match the theory better to the actual realistic situa- 
tion, we performed simulations with multiple atom pairs. 
We have two elongated volumes, with the atoms sparsely 
distributed, giving varying distances for the atom pairs. 
V is also angular dependent, this complicates the prob- 
lem a lot, also because jJsp and fidp have several possible 
polarizations. We will assume the angular dependence is 
averaged out and we use for Vdd 



(a;2 + {y + dy + z^) 



3/2 ' 



(55) 



with d the distance between the two cylinders. 

We randomly placed one particle in an elongated ellip- 
soid, Gaussian in 3D with full length 400 ^m and radius 
8 /im (both numbers are 1 / -^/e widths) at 25 /xm distance 
of a probe particle. For the interaction strength we used 
Eq. [34] with simply equal to Vdd (Eq. |55|. For 

MsdMpp used SOO^Ope^. We varied the mixing angle, 
while keeping the effective field (Eq. |25]) constant on res- 
onance for different n and aj=8 MHz. The mixing angle 
is defined as 



arctan 



/ Fri 



\V2Fs 



(56) 



In Fig. runs from to 7r/2, where corresponds to 
the horizontal axis and tt/2 corresponds to the vertical 
axis. We evaluated the |pp) probability after 20 /is for 
10000 random cases. 
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The result is shown in Fig. 
the red dashed lines, where only the height is fitted to 
the experimental data. 



D. Results 

As a first experiment we measured the 49p fraction, 
which is a measure for the dipole-dipole interaction 
strength, as a function of static field, for various values 
of the RF field, given in Fig. [T4j The measured data are 
depicted in black and are shown on top of a diagram of 
the resonances, like in Fig. [8] The red dashed lines be- 
long to the transition in eq. |52] and the blue dash-dotted 
lines belong to eq. [53] The green dotted lines depict the 
exact RF amplitudes of the performed scans, which are 
Frf = 0.049, 0.17, 0.4 and 0.585 V/cm. The green dots 
indicate where the resonances are expected, i.e. where the 
green dotted lines cross the red dashed and blue dash- 
dotted hues. 

We observe that the positions of the peaks correspond 
nicely to the expected positions, the green dots. The ex- 
act field values F^p and Fg have been slightly adjusted 
by a multiplication factor (maximum 2% different from 
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FIG. 14: (Color online) The 49p fraction in arbitrary units, 
depicted in black, is measured as a function of static field 
for different RF amplitudes. The red dashed and blue dash- 
dotted resonance lines correspond to respectively eq. |52| and 
53 and run from the -1-3 photon transition (lower left corner) 
to the -14 or -11 photon transition (upper right corner). The 
0-photon transition is depicted with thicker lines. The reso- 
nance peaks are expected where the green dotted lines cross 
with the red dashed and blue dash-dotted lines, indicated with 
green dots. 



1) to fit the expected positions, because the exact effec- 
tive plate distance is not accurately known. The mea- 
surements show that indeed more sidebands occur for a 
stronger RF field, or in other words, more multi-photon 
transitions can be observed. With an amplitude of only 
0.049 V/cm just two sidebands are populated, the one 
above and the one below the original states. We can also 
say that we have observed the n — l transition, where the 
two-atom system absorbs one RF photon, and we have 
the n = transition, the original resonance, and there is 
the n = — 1 transition, where the two-atom system emits 
one RF photon through stimulated emission. When we 
increase the RF amplitude to 0.17 V/cm we observe side- 
bands from n = 2 down to n = — 5. Surprisingly, at this 
field the original state is completely depleted: the n = 
peaks are missing. At Frf — 0.4 V/cm the peaks around 
n = —4: are missing, while at Frf = 0.585 V/cm the peaks 
around n — —9 have disappeared. All these missing peaks 
are nicely predicted in Fig. [8] by the green, red and or- 
ange areas or the white parts in the resonance curves. 
Here the generalized Bessel function and therefore the 
strength of the transition goes to zero, or in other words 
destructive Stiickelberg interference occurs. 

As a second experiment, we measured the 49p fraction 
as a function of RF amplitude Fjip for ^5=0. In other 
words, we made a vertical cut through Fig. |8] along the 
vertical axis. The results are depicted in Fig. [l5] The 
observed peaks fit the expected resonances, depicted as 
red dashed lines (transition in Eq. 52 1 and blue dotted 



lines (transition in Eq. 53 ) well; again a small adjustment 
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FIG. 15: (Color online) The 49p fraction as a function of RF 
amplitude at zero static field. The red dashed and blue dotted 
vertical lines depict the expected positions of the even photon 
transitions for both transitions given in Eq. [52] and |53| 



of the field has been made. This figure shows in fact the 
pure AC-Stark shift of the states, which is \ceF^p (see 
Eq. 24 1 . Interestingly, we observe only even-photon tran- 



sitions. In the absence of a static field only an even num- 
ber of photons can be transferred, due to the selection 
rules of photon transitions, as was the case in [23^ and is 
also visible in Fig.|8] Both atoms make a A£ — ±l transi- 
tion, and one photon - with angular momentum 1 - goes 
from one atom to the other. Adding an even number 
of photons the transition is possible, because then the 
total angular momentum of these RF photons can add 
up to zero. With an odd number of RF photons, how- 
ever, some angular momentum remains, and the photons 
can not be absorbed by the two atoms. Already in a 
small static field odd-photon transitions are allowed, as 
can be seen in Fig. |8] because the static field can con- 
tain some angular momentum (but no energy). In other 
words, with a static field other £ states are mixed in the 
atomic wavefunctions, as described in section IIII A| and 



the photons can couple these parts of the wavefunction. 
Note that only a tiny fraction of other £ states is already 
enough to break up the selection rules. In other words, 
for a small field the radial part of the wavefunction is 
hardly changed, but the time dependent part is changed 
significantly. 

A third experiment illustrates the Stiickelberg oscil- 
lations most clearly. Here we keep the effective field 



(Eq. 25 1 constant and fixed to a resonance and we vary 
the mixing angle (Eq. [56]). This requires adjustments of 



both the field amplitude and the static field. In Fig. [T6| 
the result is depicted for mixing angles between and 
tt/2 for the n = — 3 to ti = +2 photon transitions of 
Eq. |52[ In fact the measurements run over the gray lines 
in Fig. |8] starting on the bottom. The clear presence 
of the Stiickelberg interference patterns in the data, es- 
pecially at the positions of the destructive interference. 




n/8 nl4 3nl8 

arctan[FRF/(V2Fs)] 



FIG. 16; (Color online) The 49p fraction measured for dif- 
ferent mixing angles (Eq. 56 1. From left to right we have 



increasing Frf and decreasing Fs- The effective field is kept 
constant at the n- photon resonance for each plot. The red 
dashed lines show the result of the simulation discussed in 
section [III C| The blue dash-dotted lines depict the squared 
generalized Bessel function (Eq. 23 1 directly, which can be 
seen as a small transfer approximation. The green dotted 
vertical lines depict the boundaries of the classically allowed 
region and correspond to the lower blue dashed lines in Fig. [8] 



proves that our system is coherent for at least one period 
of the RF oscillation, i.e. 125 ns. Since the crossing is tra- 
versed mostly diabatically, many oscillations are needed 
to build up the pp population. It can not be proven from 
these data if this build up occurs coherently or incoher- 
ently. However, if we look at Figs. [14] and [15] we can 
estimate the ratio of peak width to peak distance to be 



approximately 5. This means that N in Eq. 49 is at least 
5, so we can say that the interaction is coherent for at 
least 5 RF periods, i.e. 0.6 fis. In reality the coherence 
time could be longer, because the resonance peaks are 
probably broadened through other mechanisms. 

The red dashed lines in Fig. 



16 show the results of 



the simulations, described in section jlll C] Here only the 
height and the vertical offset are adjusted to fit the data. 
The simulations fit the data nicely, proving that this 
straightforward calculation reflects the real experimen- 
tal situation quite well. The blue dash-dotted lines in 
Fig. [TH] depict the squared generalized Bessel function di- 
rectly (Eq. 23 1 . This reflects a situation of weak or short 
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interaction, since the population of the pp-state then be- 
comes Ppp « V^J^t^. When the interaction is stronger, 
the pp probabihty oscillates between and 1: the pro- 
cess is saturated. However, because multiple distances 
occur and therefore different interaction strengths, these 
oscillations are washed out. The result is that the peaks 
of the blue dash-dotted lines are broadened, as visible in 
the red dashed lines of the simulation. We see again that 
for small angles (DC-dominated cases) the population is 
fully in the n = state. With increasing RF amplitude 
the population in this state decreases and more and more 
higher order states are populated. For higher \n\ the start 
of the increasing population occurs at a larger RF am- 
plitude (or larger mixing angle). This is also visible in 
the classical boundaries, depicted as green dotted lines, 
equivalent to the blue dashed lines in Fig. [8] Beyond 
these boundaries we observe the Stiickelberg oscillations. 
On the right hand side of the plot we see again that the 
odd photon transitions go to a minimum and the even 
photon transitions go to a maximum, because on this 
side the static field is zero and the selection rules apply 
as explained in connection with Fig. |15| 

IV. CONCLUSIONS 

We have studied both experimentally and theoretically 
interactions between atoms under the influence of an os- 
cillating electromagnetic field. States with a quadratic 
Stark shift (polarizable states) show a fundamentally 
different behavior than states with a linear Stark shift 
(states with a permanent dipole moment). A compari- 
son between the two cases has been made using the Flo- 
quet approach. Where in the linear case the resonance 
frequency depends purely on the static field and the cou- 
pling strength depends purely on the amplitude of the 



field, in the case of the quadratic shift, both the res- 
onance frequency and the coupling strength depend on 
both the amplitude and the static offset of the field, and 
no separation of variables can be made. For the coupling 
strength in the quadratic case we invoked the generalized 
Bessel function. Its role is similar to the regular Bessel 
function for the linear case. The behavior of both func- 
tions is studied as a function of frequency and compared 
to the classical limit; the energy values that occur the 
most in the oscillating field correspond to the resonances 
with the strongest coupling strength, apart from some os- 
cillatory behavior. The oscillatory behavior of the Bessel 
function as well as the generalized Bessel function can 
be explained in terms of Stiickelberg oscillations, which 
is an interference effect between the developed phases of 
the two interacting states. 

In the experiment the resonance positions and the in- 
teraction strength fully fit the described theory. Destruc- 
tive Stiickelberg interference is clearly observed, which 
proves that our system is coherent for at least one period 
of the RF oscillation, 125 ns. From the ratio of the mu- 
tual distance the peaks and the width of the peaks, the 
lower bound is increased by a factor of 5 to 0.6 /xs. Sim- 
ulations with multiple atom pairs at different distances, 
corresponding to our experimental situation of two cylin- 
ders of atoms, fit the observed interference patterns well. 
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